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I— > ' 

I provide a brief summary of the theory presentations at ISMD 2011 related to the 
phenomenology of small- a? QCD evolution and its application to particle production and 
fluctuations in high-energy hadron and heavy-ion collisions. I also mention some challenges 
for quantitative phenomenology which emerged from the LHC, such as understanding the 
long-range "ridge" in high-multiplicity p + p collisions, the transverse momentum distribu- 
O |' tions in p + p at semi-hard p± , and the origin and scale of density fluctuations in the initial 

I ' state of heavy-ion collisions. 

(— i ' §1. Introduction 

| The physics of multi-particle dynamics has a long history. Indeed, ISMD 2011 

on Miyajima-Island, Hiroshima, was the 41 st meeting in the International Sympo- 
sium on Multiparticle Dynamics series. The talks have mostly focused on physics 
at accelerator facilities, in particular at the BNL-RHIC and CERN-LHC colliders. 
They covered a very broad range of topics including particle production and multi- 
. particle correlations in high-energy collisions, studies of hard QCD jets and of the 

underlying event, collective effects in QCD at high temperatures and baryon den- 
sities, femtoscopic Bose-Einstein correlations, as well as electroweak physics of the 
standard model and beyond. However, we had talks also on cosmic ray physics and 
its relation to accelerator based high-energy QCD, and finally on strong-field QED 
^ ■ with high intensity lasers. It is, indeed, a time of exciting progress in this field and 

ISMD 2011 provided a nice overview as well as a platform for stimulating discussions. 

This write-up is intended to offer a brief summary of the theory presentations 
at ISMD 2011 related to small- £ QCD; more specifically, of the phenomenology of 
small- x QCD evolution and its application to particle production and fluctuations 
in high-energy hadron and heavy-ion collisions. 

§2. Phenomenology of running-coupling BK evolution 

Collisions of protons or heavy ions such as lead at high energies release a very 
large number of gluons from their wave functions which dominate the production of 
new particles and anti-particles. In fact, the wave function of a hadron boosted to 



*' e-mail: Adrian.Dumitru@baruch.cuny.edu 



typeset using VTpT^.ds (Ver.0.9) 



2 



A. Dumitru 



(nearly) the light cone is so densely packed with gluons that they may "overlap", 
leading to non-linear interactions. 1 ) Therefore, at high energies the colliding hadrons 
can be treated as a high occupancy gluon field. This dense system is nowadays 
referred to as a Color Glass Condensate (CGC). For a recent in-depth review of the 
CGC we refer to ref. 2 ) 
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Fig. 1. Left: Albacete's unintegrated gluon distribution &(k±,x) obtained from rcBK evolution of 
the MV-model dipole, eq. (2-4). Right: associated yield in CMS ridge peaks about a semi-hard 
scale. 



The evolution of observables with energy is described by "small-x" QCD evo- 
lution equations (where x ~ l/V 5 ^ s inversely proportional to the energy of the 
beams). The two-point function of two light-like Wilson lines, in particular, 
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(2-1) 



determines DIS structure functions and single-inclusive cross sections in hadronic 
collisions*^ . Af(r, x) represents the scattering amplitude of a color dipole of transverse 
size r. Its evolution with energy is governed by the non- linear so-called Balitsky- 
Kovchegov (BK) equation, 

dAf(r x) f 

— — ; ' / = / d 2 n K(r,r 1 ,r 2 ) [Af(ri,x) +Af(r 2 ,x) -J\f(r,x) - N{n,x) N{r 2 ,x)} 
dlog(x /x) J 

(2-2) 

Present numerical solutions (and corresponding applications) of this evolution equa- 
tion employ a kernel valid to running coupling accuracy, 4 ) 
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(2-3) 



From the Fourier transform of the scattering amplitude for a dipole in the adjoint rep- 
resentation, J\fA(r,x) = 2N(r,x) — N 2 (r,x), one obtains the "(dipole) unintegrated 



On the other hand, the connected contribution to multi-particle production in general involves 
higher n-point functions. 
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gluon distribution" <P(k±,x) ~ k\ N(k±,x) which determines the single-inclusive 
gluon production cross section. Owing to the non-linear term in (2-2) this function 
vanishes as k± — > and peaks about a semi-hard scale ~ Q s (x), as shown in fig. 1. 
Although the direct analogy (so far) is somewhat murky, it is nevertheless tempting 
to show it side by side with the "ridge yield" obtained by CMS for high-multiplicity 
p+p collisions at 7 TeV 5 ^ (also shown by J.-H. Kim at this meeting); the associ- 
ated yield in the ridge evidently peaks about a semi-hard scale, too, which needs to 
be understood more clearly (for some initial ideas to relate the CMS ridge to the 
properties of non-linear evolution see refs. 6 )). 

As shown in the talks by Albacete, Fujii and Nara, the rcBK solution also 
provides a good description of HERA DIS structure functions at low x, of hadron 
transverse momentum distributions in the forward region of p+p and d+Au collisions 
at RHIC, and of pj_-integrated hadron multiplicities in Pb+Pb collisions at the LHC. 
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Fig. 2. Left: AAMQS rcBK fit to HERA structure functions. 7 ^ Right: charged hadron multiplicities 
versus centrality for Pb + Pb collisions at 2.76 TeV. 8 ' Presented by J. Albacete. 



Some examples are shown in figs. 2 and 3. It is worth noting though that the AAMQS 
DIS fit requires an initial condition for eq. (2-2) which is not like the "MV model" 
dipole, 

A/W(r, x ) = 1 - exp (~ r 2 Q 2 s {x ) log , (2-4) 

but one which falls off more rapidly at small r: 

MMvAr,xo) = l-exp^-^(r 2 Q 2 s (x )y\og^j , ( 7 > 1) . (2-5) 

This finding is underlined by the charged hadron transverse momentum distributions 
in p + p collisions measured at the LHC, see fig. 4. The transverse momentum 
distributions at p± > Q s have been obtained by first computing the transverse 
momentum distribution of gluons from the fc^-factorization formula which is then 
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Fig. 3. Hadron transverse momentum distributions in the forward region of p + p and d + Au 
collisions at RHIC as presented by H. Fujii at this meeting. 9 ' 
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Fig. 4. Hadron transverse momentum distributions in the central region of p + p collisions at 
7 TeV. 10) Theory K factor adjusted to fit data at px = 1 GeV. CMS data from. 11 ' 



convoluted with a DSS gluon — > hadron QCD fragmentation function: 

dN h ± f dz dN g 



.12) 



Here, pj_ = Z Q± is the transverse momentum of the produced hadron and the scale 
Q 2 may be chosen as q\. For the first time, the extreme energy reach of LHC enables 
us to perform quantitative tests of the theory at high transverse momentum p± but 
small x. Evidently, the MV model initial condition with rcBK evolution (plus k± 
factorization) disagrees with the CMS data. On the other hand, it is comforting to 
see that the AAMQS-modified MV 7 initial condition which worked for F% at HERA 
provides a good match to the LHC high-pj_, low- cc data as well. 

As is often the case, answers lead to more questions. Such as: is there any theory 
which would at least motivate the successful AAMQS guess (2-5) ? The MV model 
initial condition (2-4) can be derived from an effective theory treating the large- x 
valence charges as recoilless sources on the light cone (analogous to the Weizsacker- 
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Williams equivalent photon approximation); in the limit of a very large number of 
such valence charges, color fluctuations are described by a quadratic action 13 ) 

d 2 x ± ^ p a p a , (2-7) 

where p 2 (x±) ~ N va i is proportional to the number of valence charges at a given 
transverse coordinate. For a nucleus, N va i ~ A^A 1 / 3 is proportional to the number 
of colors and to the longitudinal thickness. 

The MV 7 initial condition (2-5), on the other hand, has not been derived (so 
far) from an effective action for the large- a; valence charges. Worse yet, there is no 
theoretical understanding of the dependence of the parameter 7 on N va i resp. on the 
thickness A 1 / 3 of a nucleus; 7 would have to be re-fitted for every single hadron or 
nucleus. Understanding the A dependence of the AAMQS parameter 7 will also be 
important for 

i) the ratio of pA to pp spectra at the LHC: some estimates using rcBK unintegrated 
gluon densities suggest that R p a might be suppressed at LHC energies, even away 
from the proton fragmentation region, due to the effects from small-x evolution; 15 ) 
and 

ii) for predictions or fits of structure functions of heavy ions at small x which will 
hopefully be measured by a future electron-ion collider elC. 14 ) 

One possibility is that the AAMQS initial condition arises from corrections to 
the MV model action which involve higher powers of the color charge density p. 
Such terms are allowed by SU(N) invariance and involve couplings suppressed by 
additional powers of 1/gA 1 / 3 : 



= I d 2 x ± \^p a p a - -d abc p a p b p c + -p a p a p b p b 
J yip, 1 k 3 re 4 



(2-8) 



Here, p 2 ~ g 2 A 1 / 3 , K3 ~ g 3 A 2 / 3 , and K3 ~ g 4 A 16 ) This action would allow for 
definite predictions for the A-dependence of the initial condition for rcBK evolution. 
It will be interesting to see if (2-8) is able to reproduce the essential features of the 
AAMQS MV 7 model. 
Coming back to 

dN pA /d 2 p± 

RpA ^ = iVcon dN pp /a- P± ' (2 ' 9) 

I. Potashnikova and B. Kopeliovich in their contribution pointed out that there are 
some subtleties in defining the "number of binary parton-parton collisions" N co \\. 17 ^ 
This is commonly taken as 

pp Irp \ 

N co n = ^£f, (2-10) 

where T4 denotes the thickness function of the target nucleus, Pi n is the probability 
for an inelastic pA interaction, and (•) denotes an average over events. The authors 
argued that the cross section in eq. (2T0) should be corrected for the fraction of 
missed diffractive events. If this represents a significant fraction then N co \\ is reduced 
and this would affect our interpretation of nuclear effects as measured by R p a- 
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Fig. 5. Cross section for two jets with p± > 10 GeV separated by a rapidity interval Ar\ and 
azimuthal angle A<j>\ p + p collisions at y/s = 7 TeV. From ref. 18 ' 

Another important test for our present knowledge of the unintegrated gluon 
densities has been emphasized by K. Kutak at this meeting. That is, the single- 
inclusive cross section at non-central rapidities where the evolution is asymmetric; 
also, the angular distribution of semi-hard jets separated by a rapidity interval Z\ry 18 ) 
shown in fig. 5. 

§3. Density fluctuations in p + p and A + A collisions 

There has been a lot of interest recently in understanding the magnitude and 
length scale of fluctuations of the density of produced gluons in the transverse plane; 
several presentations at this meeting addressed related questions. In their talks, 
A. Mocsy and P. Sorensen noted the analogy to the power spectrum of CMB fluctu- 
ations which provides information on the evolution of the early universe. An artists 
view is shown in fig. 6. They presented the harmonic spectrum of the initial density 
distribution (preceding the hydrodynamic expansion in A + A collisions) in terms of 
the eccentricities 

2 (r 2 cos n(f)) 2 + (r 2 sinm?i>) 2 

(•} denotes an average over the initial distribution of produced gluons in the trans- 
verse plane, 20 ) dN/drjd 2 r. One commonly considered source of inhomogeneities is 
due to the fluctuations of the positions of the nucleons in the nuclei before the colli- 
sion. In the absence of a length scale (such as a non-zero radius of the nucleon) for 
these fluctuations, and for central collisions, Sorensen argues that the spectrum of 
e n is flat*) . 



The peak of e| in n g- 6 for b = collisions is due to the intrinsic deformation of a Au nucleus. 
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Fig. 6. Left: Artists view of the expansion of the universe and of the survivial of initial quantum 
fluctuations as CMB fluctuations measured today; and of the expansion of fluctuations in p + p 
and A+A collisions. 19 ' Figure by A. Doig. Right: harmonic spectrum of density inhomogeneities 
in central A + A collisions for point-like nucleons. Presented at ISMD 2011 by A. Mocsy and 
P. Sorensen. 
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Fig. 7. Right: Eccentricity spectrum averaged over Monte-Carlo ensembles of nucleon configura- 
tions with Gaussian valance parton distributions, obtained from k± factorization with rcBK 
unintegrated gluon densities. Presented by Y. Nara. 21 ' 



Y. Nara showed the Fourier spectrum of eccentricities obtained from a simulation 
of heavy-ion collisions based on the k± factorization formula with rcBK unintegrated 
gluon densities, averaged over an ensemble of configurations of the colliding nuclei. 
Again, here the only source of fluctuations are the random positions of the nucleons 
before the collision (the average is a Woods-Saxon distribution); they determine the 
locations of the valence charges for small-x evolution. The valence charge profile 
function was assumed to be Gaussian rather than a delta- function. The Fourier 
spectrum of eccentricities then drops towards higher harmonics, the faster the bigger 
the Gaussian width B. A careful understanding of higher harmonics, in particular 
of those moments which are dominated by fluctuations, will require constraints on 
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the Gaussian width B from exclusive vector meson production in DIS, and from 
multiplicity distributions in p + p and p + A; see, e.g., ref. 22 ) and references therein. 

Fluctuations of the nucleon configurations in the colliding nuclei are very likely 
not the full story, however, as has been pointed out by C. Flensburg and Y. Hatta. 
There are intrinsic fluctuations of the number of emitted gluons even for a fixed 
number of longitudinally stacked projectile and target nucleons (or "participants") 
which are actually amplified in a collision of strong chromo fields. 23 - 1 Furthermore, 
one expects fluctuations due to stochastic dipole splitting events as implemented in 
Flensburg's DIPSY Monte-Carlo; 24 ) the probability distribution for splitting is given 
by 

N c (x-y) 2 
' 2ir 2 (x — z) 2 (z — y) 5 

if a fixed coupling constant a s is assumed. Hatta stressed that (3-2) results in large 
fluctuations of the number of dipoles which are strongly correlated in the impact 
parameter plane. 

In contrast, so far solutions of the rcBK equation have been averaged over all 
random emissions, and corresponding particle production models have not considered 
the effects from this type of fluctuations. Fig. 8 shows that fluctuations in dipole 
splittings can increase £3 quite significantly, while €2 appears to be largely unaffected. 



dP = as7 ^ J ^A_y> i SzdY , (3-2) 
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Fig. 8. Left: e n from the Monte-Carlo stochastic dipole splitting model DIPSY, from a geometric 
"wounded nucleon / Glauber" model, and from Nara's MC-KLN model 26 ' (fluctuations of large- 
x sources only). Presented by C. Flensburg. Right: DIPSY prediction for the fourth cumulant 
of non-flow correlations for high- multiplicity p + p collisions at 7 TeV. Preliminary ALICE data 
from. 27 ' Presented by Y. Hatta. 



The presence of such large fluctuations in the gluon density is evident from the 
multiplicity distribution in p+p collisions at high energies. It also makes it feasible to 
define Fourier "eccentricities" of the transverse area occupied by produced gluons in 
high- multiplicity p+p collisions, in analogy to heavy- ion collisions. What is measured 
though are the momentum- space azimuthal anisotropics v n in the final state. DIPSY 
predicts that the large fluctuations would lead to ffl^l (measured via four-particle 
cumulants) which is negative, corresponding to "non-flow" correlations. 25 ^ 
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§4. Summary 

With the ongoing program (and possible upgrades) at RHIC and with the recent 
advent of LHC, this is clearly an excellent time for progress in high-energy QCD. 
Tools have been constructed, and constantly refined, to address questions about the 
strong color fields of hadrons and nuclei boosted nearly to the light cone and their 
manifestation in particle production and in fluctuations analogous to those at the 
Big Bang. Hopefully measurements at non-central rapidities as well as for p + A 
collisions will be possible at the LHC to probe QCD in a highly non-linear regime. 
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